Squeezed Correlations and Spectra for Mass-Shifted Bosons 
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Huge back-to-back correlations are shown to arise for thermal ensembles of bosonic states with 
medium-modified masses. The effect is experimentally observable in high energy heavy ion collisions. 
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Introduction — The Hanbury Brown - Twiss (HBT) 
or Goldhaber - Goldhaber - Lee - Pais (GGLP) effect 
jl],D is widely used in heavy ion physics to measure the 
space-time geometry of such reactions. The enhanced 
correlations of bosons in outgoing states with small rel- 
ative momentum provide a Fourier transformed picture 
of the system at freeze-out. The scales measured by the 
HBT effect coincide with the lengths of homogeneity jj) . 

In this Letter we consider the effect of possible mass 
shifts in the dense medium on two boson correlations 
in general. Thus far medium modifications of hadron 
masses have been mainly considered in terms of ef- 
fects on such observables as dilepton yields and spectra. 
Hadron mass shifts are caused by interactions in a dense 
medium and therefore vanish on the freeze-out surface. 
Thus, a naive first expectation is that in-medium hadron- 
modifications may have little or no effect on two boson 
correlations, and so the usual HBT effect has been ex- 
pected to be only concerned with the geometry and mat- 
ter flow gradients on the freeze-out surface. However, in 
this Letter we show that an interesting quantum mechan- 
ical correlation is induced due to the fact that medium 
modified bosons can be represented in terms of two-mode 
squeezed states of the asymptotic bosons, which are ob- 
servables. As a by-product, we solve the finite-size prob- 
lem for two-mode squeezing, by formulating the theory 
of intensity interferometry for a finite, inhomogeneous, 
squeezed and expanding medium, which is a new result 
also for quantum optics H . 

In this Letter we assume the validity of relativistic hy- 
drodynamics up to freeze-out. The local temperature 
T(x), chemical potential pb{x), and flow fields it M (x), are 
given, for example, by those taken from ref. ||. In rel- 
ativistic heavy ion collisions, it has been observed that 
the one particle spectra can be described by thermal dis- 
tribution fairly well ||. We assume that the sudden 
(non-adiabatic) approximation is a valid abstraction in 
describing the freeze-out process in relativistic heavy ion 
collisions quantum mechanically |Q and that there ex- 
ists an abrupt freeze-out surface, E M (a;). However, the 
effects of fluctuations of that freeze-out process will also 
be considered below. 

Consider the following model Hamiltonian, 



Ho 



H = H 



dxf^ + lV^r+mg^), 



(2) 



where H is the asymptotic Hamiltonian, in the rest 
frame of matter. The scalar field </>(x) in the Hamil- 
tonian H corresponds to quasi - particles that propagate 
with a momentum-dependent medium-modified effective 
mass, which is related to the vacuum mass, mo, via 

ml(\k\) = m 2 -6M 2 (\k\). 

The mass-shift is assumed to be limited to long wave- 
length collective modes : 5M 2 (|k|) < to 2 , if |k| > A s . 

We are interested in the invariant single-particle and 
two-particle momentum distributions: 



, , d 3 N . t 

Ni (ki ) = w kl -j^- = uj kl (a ki a ki ) 



(3) 
(4) 



7V 2 (ki,k 2 ) = WkjWk^a^o^a^a^), 

+ ( a L a k 2 > ( a L a kl > + ( a L a l 2 > (« k2 « kl ) > ( 5 ) 

where a k is the annihilation operator for the asymp- 
totic quantum with four- momentum = (w k , k), w k = 
to,q + k 2 (w k > 0) and the expectation value of an opera- 
tor O is given with the density matrix p as (O) = TV pO. 
Eq. (||) has been derived as a generalization of Wick's the- 
orem for locally equilibrated (chaotic) systems in refs. 



In order to simplify notation, we introduce the chaotic 
and squeezed amplitudes, defined, respectively, as 



G c (l, 2) = V^ kl w k2 (a ki a k2 
G s (l,2) = v/WkiWka (a ki a ka 



(6) 
(7) 



Usually, the chaotic amplitude, G c (l,2) = G(l,2) is 
dominant, and carries the Bose-Einstcin correlations, 
while the squeezed amplitude, G s (l,2) vanishes: 



G 2 (k!,k 2 ) 



jV 2 (ki,k 2 ) 
JVi(ki)iVi(k 2 ) 



= 1 



|G(1,2)| 2 
G(1,1)G(2,2)' 



(8) 



dxdy0(x),5M 2 (x-y)0(y), 



(1) 



The exact value of the intercept, G 2 (k, k) = 2, is a char- 
acteristic signature of a chaotic Bose gas without dy- 
namical 2-body correlations outside the domain of Bose- 
Einstein condensation fnO|. 
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For a hydrodynamic ensemble, eq. (g) reduces to the 
special form derived by Makhlin and Sinyukov 



V 



G(l,2) 



J-3 J d^K^e^ * n(x,K u 



(27T) 



(9) 



Note that the relative and average pair momentum co- 
ordinates, q. 2 = wi — u>2, qi,2 = ki — k2, and Ki :2 = 
|(ki + k 2 ) appear in (||). These variables arise natu- 
rally whenever the Wigner operator is used, even in to- 
tally non-equilibrium semi-classical limits |n[]. The va- 
lidity of the approximations leading to (|9j) requires the 
width of (7(1, 2) as a function of the relative momentum, 
q = |qi,2 1 , to be small. That width is given by ~ 1/i?, 
where R is a characteristic dimension of the system. The 
semi-classical limit corresponds to KR ^S> 1, where if 
is |Ki f 2|. Note that y/oJj^ cuk 2 ~ ^12 m this case. For 
gi? < 1, the second term in (^J) describes the minimal 
quantum interference associated with the indistinguisha- 
bility of the bosons. The integration over the freeze-out 
surface, ^(x), is implemented with the invariant mea- 
sure d 3 'E fl K^ 2 that reduces to K°d 3 x in the special case 
of a constant freeze-out time. 

Results for a homogeneous system - The terms ne- 
glected in (||) involving G s (l,2) become non-negligible 
when mass shift becomes non- vanishing, i.e., <5M 2 (|k|) ^ 
0. Given such a mass shift, the dispersion relation is mod- 
ified to il k — oj^ — SM 2 (\k\), where is the frequency 
of the in-medium mode with momentum k. The anni- 
hilation operator for the in-medium quasi-particle with 
momentum k, &k ; and that of the asymptotic field, ak, 
are related by a Bogoliubov transformation [12j : 



s b 1 



Ci+sii, (10) 

where Ck = coshfrk], Sk = sinh[rk] and rk is given by 
1 



rk 



■ log(w k /fik) 



(11) 



We introduce the shorthand, C x and S_ x , to simplify 
later notation. As is well-known, the Bogoliubov trans- 
formation is equivalent to a squeezing operation, and 
so we call rk the mode dependent squeezing parameter. 
While it is the a-quanta that are observed, it is the b- 
quanta that are thermalized in medium. Thus, we con- 
sider the thermal average for a globally thermalized gas 
of the 6-quanta, that is homogeneous in volume V: 



f ' = Z ,,X!, l4(2^ 



When this thermal average is applied, 

G c (l,2) = V^kl [(C\C 2 ) + {S_.SU) 
G s (l,2) - V^kT \(SUC 2 ) + (CjSla) 



(12) 



(13) 
(14) 



If this thermal b gas freezes out suddenly at some time 
at temperature T, the observed single a-particle distri- 
bution takes the following form: 



Wk^i(k), 



m(k) = |c k | 2 n k + |s_ k | 2 (n-k + l), 
1 



exp(r» k /r) - 1 



(15) 

(16) 
(17) 



This spectrum includes a squeezed vacuum contribution 
in addition to the mass modified thermal spectrum. Its 
shape, however, remains essentially that of the conven- 
tional thermal distribution. Although the squeezed vac- 
uum contribution results in a slowly decaying power- 
law tail of the single-particle spectra, this power-law tail 
arises for typical values of mass-shifts and temperatures 
at high transverse momentum only, e.g. above |k| > 1 
GeV, where the local hydrodynamic description breaks 
down and correlated pQCD (mini-jets) starts to play a 
dominant role. We emphasize that the spectra and cor- 
relation functions should be described by the same set of 
model parameters, and any signal found in the correla- 
tion function could be cross-checked against its contribu- 
tion to 2Vi(k). 

In the homogeneous limiting case, G c (l,2) cx V5i^, 
while G s (l,2) cx VSi^ 2 - The resulting two particle cor- 
relation function is therefore unity except for the parallel 
and antiparallel cases: 



G 2 (k,k) = 2, 
G 2 (k,-k) = H 



Ifkfk^k + c_ k S-k( n -k + X )! 2 
7ii(k)m(-k) 



(18) 
(19) 



The dynamical correlation due to the two mode squeez- 
ing associated with mass shifts is therefore back-to-back as 
first pointed out in [^2[ . The HBT correlation intercept 
remains 2 for identical momenta. Evaluating eq. ( |l9| ) for 
T = 140 MeV, |k| = 0, 300 or 500 MeV for <f) mesons, as 
a function of the medium modified mj,, one finds back- 
to-back correlations (BBC) as big as 100 - 1000 for rea- 
sonable values of mj. 

It follows from eq. (pit) that the BBC are unbounded 
from above, 1 < G 2 (k, — k) < 00. As |k| — > 00, 
G 2 (k,-k) ~ 1 + l/|s_ k | 2 ^ 1 + l/ni(k) -> 00. Hence 
at large values of k, particle production is dominated 
by that of back-to-back correlated pairs for any non- 
vanishing value of in-medium mass-shifts. This huge en- 
hancement renders our effect measurable. This predic- 
tion has to be contrasted to the results of ref. 0, which 
predicted large values for the HBT correlations, G(k, k) 
as a consequence of an incorrect Bogoliubov transfor- 
mation. Ref. U discussed a different kind of back-to- 
back correlation, the particle - antiparticle correlation 
(PAC), which is based on a classical current formalism, 
and showed C 2 PAC \\s., — k) < 3. However, their upper 
limit was shown to be reachable for massless quanta at 
|k| = only. For massive final state bosons and for any 
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realistic finite duration of particle emission, the effects 
found in ref. || are suppressed by at least 2 -3 orders 
of magnitude, resulting in an unmeasurably small effect. 
In contrast, we consider here an entirely quantum effect, 
which cannot be obtained in classical current formalism. 
The unlimited strength of the back-to-back correlations of 
mass-shifted mesons in our case results in a measurable 
effect even if the decay of the medium is not completely 
sudden and orders of magnitude suppressions reduce its 
strength. 

Suppression by finite emission times. To describe 
a more gradual freeze-out, the probability distribution 
F(ti) of the decay times ti is introduced. (The sud- 
den approximation is recovered in the F(ti) = S(U — to) 
limiting case.) The time evolution of a k (i) will be 
flk(i) = Ok(ii) exp[— iujk(t — ti)], which leads to 



G 2 (k,-k)-l = 

= l c k s k"k + cl k a_ k (n_ k + 1)| 2 
rc.i(k)rci(-k) 



|^K+^-k)| 2 , (20) 



where F(ui) — J dtF(t) exp(~iujt). For a typical expo- 
nential decay F{t) = 9(t — to)T exp[— T{t — to)], the sup- 
pression factor is 

\F(u k + w _ k )| 2 = i/[i + K + ^_ k ) 2 /r 2 ]. 

In the adiabatic limit, r — * 0, this factor suppresses 
completely the BBC, while in the sudden approximation, 
r — > oo, the full BBC are preserved. For St = K/T = 2 
fm/c, and for BBC of <f> mesons with uj ~ I — 2 GeV, 
one finds that the BBC are suppressed by factor ~ 10~ 3 . 
As shown in Fig. 1, the BBC survive this large sup- 
pression with a measurable strength, as large as 2-3, the 
scale of HBT correlations. This emphasizes the enormous 
strength of squeezed BBC for mass-shifted bosons. 

Results for inhomogeneous systems. Following Ref. || , 
we divide the inhomogeneous fluid into independent cells 
labeled i and assume 



n& = n^ ex p 



-{Hi - jUiJVO/Ti 



(21) 



where pi is the local thermal density matrix of cell i at 
the decoupling time ti, and AT, is the local number oper- 
ator of b quanta. The hydrodynamic limit applies only if 
the cells are small compared to the scale of change of the 
temperature, flow and chemical potential fields, however, 
they are big enough so that one can apply quantization 
and statistics cell by cell. The validity of hydrodynamics 
in heavy ion collisions is not obvious, but must be checked 
in each reaction through the consistency with both single 
and double inclusive spectra. In each cell, the field can be 
expanded with creation and annihilation operators, and 
Hi is diagonalized by a local Bogoliubov transformation. 
The amplitudes @Q 

can be evaluated assuming that 
the 6-quanta satisfy locally a generalized eq. (|l0|): 



G c (l,2) = — L J d^(z)Xf ;2 e^-* 

x [|ci, 2 |V, 2 + |a_i,_2| 2 (n_i,_2 + 1)] , (22) 
G.(l,2) = — L | dV^iV^ 1 ' 2 " 

X [s* lj2 C2,-in_i )2 + Cl,_2S* 2 ,l(«l,-2 + 1)] ■ (23) 

Here d A a^{x) = d 3 S^(x;r/) F(rf)drf is the product of 
the normal-oriented volume element depending paramet- 
rically on tj (the freeze-out hypersurface parameter) and 
the invariant distribution of that parameter F{jf), and 
u^(x) is the local flow vector at freeze-out. The other 
variables are defined as follows: 



Cl k (x) = \Jml - k^kf, - 8M 2 {x, k), 
k = J-k^ku, 



k^(x)=u^(x)iJ ki Tk^x), 
k*£(x)=u»(x)n ki (x)T~kT(x), 



(24) 
(25) 

(26) 

(27) 
(28) 



n id {x) = l/[eM{KZ 3 (x)u^x)~p(x))/T(x)] - l], (29) 

r(i,j,x) - - log [(K^(x))/(K*»(x)u„(x))] , (30) 
Ci tj = cosh[r(i, j,x)], s itj = smb. [r(i,j, x)}, (31) 

where i,j = ±1,±2 and the mean and the relative mo- 
menta for the a(6)-quanta are defined as K^*j^(x) = 

[k ( *^(x) + k^(x)}/2 and q^J" = k^ - kV'», respec- 
tively. We assume of course that the mass shift and hence 
squeezing is non-vanishing over only a finite domain of 
the freeze-out hypersurface. In terms of eqs. ( p2| , p3| ), 
the following new expressions describe the particle spec- 
tra and the correlation function in the presence of local 
squeezing: 



AT 1 (k 1 ) = G c (l,l), 

C 2 (k!,k 2 ) = 1 



(32) 



|G C (1,2)| 2 



|G S (1,2)| 2 



G C (1,1)G C (2,2) G C (1,1)G C (2,2) 

Fig. 2 illustrates the novel character of BBC for two iden- 
tical bosons caused by medium mass-modifications, along 
with the familiar Bose-Einstein or HBT correlations on 
the diagonal of the (k l5 k 2 ) plane. 

P article- antiparticle pairs — As the Bogoliubov trans- 
formation always mixes particles with antiparticles, the 
above considerations hold only for particles that are their 
own antiparticles, e.g. the <j) meson and 7r°. However, 
the extension to particle - antiparticle correlations is 
straightforward. Let "+" label particles, "— " antiparti- 
cles if antiparticle is different from particle, let "0" label 
both particle and antiparticle if they are identical parti- 
cles. The non-trivial correlations from mass-modification 
for pairs of (++), (H — ) and (00) type read as follows: 
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C++(k 1 ,k 2 ) = l 
C+-(k 1 ,k 2 ) = l 
C 2 00 (k 1 ,k 2 ) = l 



|G C (1,2) 



G C (1,1)G C (2,2) : 
|G,(1,2)| 2 

G C (1,1)G C (2,2)' 
|G C (1,2)| 2 



(33) 
(34) 



|G s (l,2)p 



G C (1,1)G C (2,2) G C (1,1)G C (2,2)' 

(35) 

where we have assumed for simplicity that mass- 
modifications of particles and antiparticles are the same. 

Summary — We formulated the theory of particle cor- 
relations and spectra for bosons with in-medium mass- 
shifts, which predicts the existence of unlimited back- 
to-back correlations of (j>4>, K + K~, 7r°7r° and 7r + 7r~ 
pairs that could be searched for at CERN SPS and up- 
coming RHIC BNL heavy ion experiments |L3|. These 
correlations not only survive orders of magnitude finite 
time suppressions, but may also be utilized to determine 
the freeze-out time distribution and in-medium hadronic 
mass modifications. 
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Fig. 1. Dependence of the back-to-back correlations 
(BBC) on the medium modified <fi meson mass, for 
T = 140 MeV and fi = 0, where solid, dashed and dotted 
lines stand for |k| = 0, 300 and 500 MeV, respectively. 
The magnitude of the BBC is large in spite of the fi- 
nite time suppression factor 1/[1 + (2Stuji ! ) 2 ] ~ 0.001, for 
St = 2 fm/c. 
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Fig. 2. Schematic illustration of the new kind of cor- 
relations for mass shifted ir° pairs, assuming T = 140 
MeV, G c - exp[-q 2 2 R 2 G /2], G s ~ exp[-2K 2 2 R 2 G ], with 
Rc = 2 fm. The fall of the BBC for increasing values of 
|k| is controlled here by a momentum-dependent effective 
mass, ml = m^l + exp(-k 2 /A 2 )] with A s = 325 MeV 
in the sudden approximation. Without the A s cutoff, the 
BBC would increase indefinitely as |k| — > oo. 
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